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Abstract 

The Ohmic conductance and current through two quantum dots in series is 
investigated for the case of incoherent tunnelhng. A generalised master equa- 
tion is employed to include the discrete nature of the energy levels. Regions of 
negative differential conductance can occur in the I-V characteristics. Trans- 
port is dominated by matching energy levels, even when they do not occur at 

the charge degeneracy points. 
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Due to the improvement of lithographical techniques on a nanometre scale in recent years 
it has become possible to study systems that were inaccessible to experimentation before. 
The most easily controlled mesoscopic systems are defined in the 2-dimensional electron gas 
of a semiconductor. By the application of gate electrodesBJl it is possible to confine the 
electron gas effectively to produce quantum wires and dots. This allows charge quantisation 
to be observed in the form of the Coulomb blockadeB^. Moreover, the importance of size 
quantisation has been shown by Reed et al. who discerned discrete states in small quantum 
dotJ. 

The effect of the charge quantisation has been studied both for a single doti and a double 
dotii. In both cases the level spectrum could be considered continuous. The effect of level 
quantisation on the Ohmic conductance through a single dot has been studied in the limit 
of weak coupling to the reservoirs. Coherentil and incoherent methods0 lead to the same 
result in this limit. 

This paper will focus on the transport properties of a double dot with discrete equally 
spaced energy levels. The limit of completely incoherent tunnelling is studied, where the 
phase-breaking rate is large and the tunnelling process is sequential. This justifies the use of 
a semi-classical method like the master equation. It must be noted that the phase-breaking 
time is still large compared to the time it takes for an electron to traverse the dot. This 
ensures that the energy levels are quantised, due to the coherence of the wavefunctions inside 
the dot. 

I. GENERALISED MASTER EQUATION 

When it is assumed that the tunnelling rates are small compared to the Coulomb energy 
and the average level spacing, then use of the master equation is justified^. In order to take 
account of the discrete nature of the energy- levels in the dot, the master equation method 
for dots with a continuous energy level spectrumi has to be generalised. The tunnelling 
rates not only depend on the number of electrons in the dot, but also on their configuration, 
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i.e. how the electrons are distributed over the available energy levels. When there is a high 
relaxation rate, then the electrons will revert to their local equilibrium distribution between 
tunnelling events. This is likely to be the case for small tunnelling barriers which cause the 
electrons to spend a long time in the dot. Since this distribution will depend only on the 
temperature, the state of the dot may yet again be described simply in terms of the number 
of electrons present in the dot. 

Let P(p, N) be the probabihty that the system is in a state which is characterised by N 
electrons in the dot which are distributed over the energy levels according to a configuration 
p. Define T^^ (N) as the tunnelling rate coefficient corresponding to the transition {p, N + 
1} — > {p', N} by an electron tunneUing through the left barrier. Tl^{N) refers to the rate 
that corresponds to the reverse process. The tunnelling rates through the right barrier are 
defined similarly. Finally, there are the relaxation rates T^p'(N) which account for the intra- 
dot transitions at a fixed electron occupation N. The rate of change of the probability of 
occurrence of the state {p, A^} is thus given by 

^pip,N)^ j:lTgriN-mp',N-i)+Tf;f{N)Pip',N+i) 

n=L,R p' 

- E T.[n'i{N)-T!g{N-l)]p{p,N) 

n=L,R p' 

+ ^ [TfP{N)P{p', N) - Tf{N)P{p, N)] (1) 
p' 

If the system consists of several dots in series, it is necessary to define the tunnelling rates 
between the dots. Define TfJ^ {Ni, Ni±i) as the tunnelling rate corresponding to the tran- 
sition {p, A^j + 1} — > {p',Ni} and {q, Ni±i} {q', Ni±i + 1} where the subscript i labels the 
site of the dot. The evolution of the multi-particle state occupation probabilities for a dot 
not neighboured by a reservoir is given by 

jP{p,N.)= E E n'&f\Nr,N,^,)P,{p\N, + l)P,^,{q,N,^,) 

Ni±i p',q,q' 

+ E E ;f (^i±i, Ni - l)P,(p', N, - l)P,±i(g, iV,±i + 1) 
-EE T^:^f{N,-l,N,^,)P,ip,N,)P,^^{q,N,^,) 

Ni±i p',q,q' 



Ni±i p',q,q' 

+ E [r/^(Ar)P,(p', TV) - T?^\N)P,{p, N)] (2) 
p' 

For both the single and the multiple dot system the steady state current can be written as 

= e E E [Tl{{N)P{p, N) - Tfj^'Pip', N)] (3) 

p,p' N 

When it is assumed that the relaxation rates are large compared to the tunnelling rates, 
then the number of electrons in the dots is the only significant variable, since knowledge of 
this quantity allows one to deduce the probabilities of the various electronic configurations. 
Hence the formalism is greatly simplified. However, the tunnelling rates between states with 
a different occupation number have to be redefined to include the weighted sum over all 
possible tunnelling paths. Define Tli{N) and Til{N) as the tunnelling rates for an electron 
tunnelling into and out of the dot through the left barrier with other electrons present in 
the dot. 

TLm = T.Tl{{N)P{p\N) (4) 

p.p' 

TiL{N) = Y.Tf({N)P{p\N + l) (5) 



p.p' 



where P{p\N) is the conditional probability that the system is in configuration p given that 
there are N electrons in the dot. Since the master equation takes only single electron tun- 
nelling events into consideration, the only contributions that need to be taken into account 
are those where p and p' differ by the occupation of one energy level. Therefore the double 
sum over p and p' can be replaced by a single sum over the single particle energy levels k. 

TLm = T.Ttm[^-mm (6) 

k 

T,:^{N) =Y.T^^{N)P{k\N) (7) 

k 

Here P{k\N) is the conditional probability that level k is occupied given that there are N 
electrons in the dot. For inter-dot tunnelling events the rates are re-expressed in terms of a 
sum over the tunnelling paths from all levels k in dot i to all levels I in dot i ± 1 



Ti,±i(iV„iV,±i) =Y,Tlt^^{N,,Ni±{)P,{k\N, + 1)[1 - P^±i{m±i)] (8) 

k,l 

For the case of negligible energy level broadening the probability Pi{k\Ni) of finding an 
electron on level k of dot i given iVj electrons can be determined from the Boltzmann 
distribution. 

P.(t|A'.)=e-''<'<".-')^^|^ (9) 

where (3 = l/ksT and ek{Ni — 1) is the energy required to add an electron at an energy level 
k when A^j — 1 electrons are already present in the dot. Zi{Ni) is the partition function for 
Ni electrons in dot i, Zf{Ni) is the conditional partition function for Ni electrons given that 
level k is unoccupied. When the conditional probabilities are substituted in equation ^ one 
obtains 

Tunnelling events between dots will in general not preserve energy, since it is unlikely that 
energy levels will line up. The high relaxation rate suggests that all tunnelling events can 
be considered inelastic. Hence 

It follows from the two preceding equations that the tunnelling rates are inelastic in the 
total free energy of the system. This is also true for tunnelling between dot and reservoir. 
This implies that the following condition holds true at zero bias for A/" dots for all possible 
occupation numbers Ni. 
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II. THE CANONICAL DISTRIBUTION 



One of the quantities that is needed to calculate the tunnelling rate between dots of specified 
occupation number is the probability P{k\N) that an energy level k is occupied given a total 
occupation of N electrons. When no relaxation of the electrons in the dot is allowed the 
occupation of the levels will be determined by the coupling to the reservoirs. However, in the 
presence of thermahsation the levels will be filled according to some equilibrium distribution 
independent of the energies at which the electrons initially entered the dot. In metals the 
energy-levels are filled according to the Fermi- Dirac distribution function. Naively one might 
expect this also to hold for discrete energy levels and this has been assumed in various papers. 
Unfortunately, this is not generally the case. In metallic systems the levels are very closely 
spaced so that the addition of an additional electron will have no effect on the distribution, 
whereas this is not true for a discrete level spectrum. In this section the correct distribution 
will be calculated for a dot with an infinite number of equally spaced energy levels, containing 
N electrons. The calculation is done for energy levels with neghgible broadening (T <^ ksT). 

Consider the partition function Z{N) for a dot containing N electrons, with energy levels 
Cfc = eo + /cAg. 



p{N} 

where the sum is taken over all realisations p{N} with N = J^kUk- Explicitly expand out 
the contribution of the lowest energy level (either rii = 1 or rii = 0) . 



When it is assumed that there is an infinite number of energy levels, then the sum over 
the occupation numbers Uk {k > 1) can be transformed to the partition function over all 
levels (including k = 1). Rewriting = e^-i -|- and relabelling Uk —>■ Uk-i, the following 
recursion relation is obtained. 



Z(N) = g-/3(»*iei+n2e2+...) 



(14) 




n2e2+n3e3+...) 



(15) 



(1-e 



'^^^^)Z(iV) = e-^'°e-^'''^^Z{N - 1) 



(16) 
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With the boundary condition given by Z{0) = 1, this is solved to give 

^ / 1 \ 

Z(iV) = e-^-on(^^^,^) (17) 

This expression for the partition function does not exphcitly include the many-body inter- 
action. This is irrelevant in this context, since the number of electrons is kept constant. 
However, when one wants to calculate the probability that the dot contains a certain num- 
ber of electrons, the Coulomb interaction obviously needs to be included. For a simple 
model of the interaction where all pairs of electrons have an associated Coulomb repul- 
sion U, the full interactive partition function is obtained by multiplication with a factor 
exp[-/?A^(iV- l)f//2]. 

Now it is possible to do a similar calculation to obtain an expression for Z^{N), the 
conditional partition function for electrons given that level k is unoccupied. However, in 
contrast to the full partition function, this requires one to solve a recursion relation in both 



parameters k and A^. This can be avoided by using equations ^ and |T0] to yield a recursion 
relation of the conditional probability P(A;|A^) directly. 

Since the full partition functions are known, this is easily solved to give 

N q 

P{k\N) = n e"^^^' (e^A.(7v+i-p) _ (^g) 

9=1 P=l 

In order to compare this with the Fermi-Dirac distribution function, one needs to consider 
a very large number of electrons in the dot. The occupation probability will depend only on 
the difference x = k — N — 1/2. 

oo 

p^(x) = ^(_i)'?-ie-^^^''("+''/2) (20) 

q=l 

For computational reasons the above equation is rewritten as 

o . ^ o ^(.^ 1/4) J ^T=o e-2/^^^(^-/2+V4)^ sinh [/?A,(g - x/2 + 1/4)] x < 
P^{x) = 2e 2 ''^> { (21) 

E^o 6-2/3^^(5+^/2+3/4)' sinh [f3A,{q + x/2 + 3/4)] x > 



The distribution function has been plotted in figure ^ for various values of In the 

metallic limit, when the level spacing is small compared to the temperature, the distribution 
tends towards the Fermi-Dirac distribution function, as expected. When the temperature 
is comparable or smaller than the level spacing, the distribution deviates significantly. Its 
limiting behaviour is well described by another Fermi-Dirac distribution function with an 
effective temperature which is half the real temperature. A slightly more accurate result for 
the occupation probability in the limit (3A^ ^ 1 is given by (note that x 7^ by definition) 



e/3Ae(a:-l/2)_,_l X < Q 



Poo(x) = <i -"^Y (22) 

g/3Ae(a: + l/2)_,_;^ 3; > 



III. TUNNELLING THROUGH A SINGLE DOT 

When the energy level broadening is negligible compared to the temperature, i.e. in the 
limit of weak coupling to the reservoirs, the density of states in the dot can be adequately 
described by a set of delta functions. The leads are assumed to be in thermal equilibrium, 
described by the Fermi-Dirac distributions and To a first approximation, the electron- 
electron interaction can satisfactorily be described using the charging model, where each 
pair of electrons has an associated Coulomb repulsion of U . Fermi's golden rule gives the 
tunnelling rates between the dot and the reservoirs for all energy levels k. 

Tt,{N) = '^h{e, + NU) (23) 
T,\iN) = ^[l-hie, + NU)] (24) 

where F^, = 27ipl\Vl\'^ is the strength of coupling to the left reservoir. The expressions for 
the tunnelling rates through the right barrier are similar. It is assumed that the quantum 
dot can be approximated by a parabolic confining potential so that the single particle energy 
levels are equally spaced by an energy A^. In figures |^ and ^ the current and its associated 
differential conductance is plotted as a function of the applied bias for a range of energy 
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level spacings. In numerical calculations it is possible to take into account more realistic 
energy level spectra, enabling a closer comparison with experiment. When it is taken into 
account that the total spin of the system can only change by 1/2 with each tunnelling event, 
it appears that negative differential conductance may occur in specific regionJiilii. 

In general a dot with closely packed energy levels yields a higher current than a dot with 
a sparse energy level spectrum, because of the higher number of current paths available. 
When Ae ^ f/ then the metallic regime is entered. 

When the energy level spacing is not negligible, the I-V characteristics are typified by 
two energy scales, the Coulomb repulsion energy U and the bare energy level spacing A^. As 
in the metallic regime, one expects the I-V characteristic to display a current step whenever 
the maximum occupation of the dot increases by one. This happens with a period U + A^, 
since an extra electron not only has to overcome the Coulomb barrier but also has to tunnel 
to the next available energy level. In addition, there is also some fine structure which has 
an associated period of A^. This is caused by the fact that an extra current path is created 
when the bias is increased by A^. 

When Tl ^ F/j the dot will be mostly maximally occupied and the most marked current 
increases occur whenever the Coulomb blockade can be overcome. This means that the 
period f/ + A^ is accentuated. In the opposite regime Tl <^Tr the dominant period is the 
level spacing. Under normal operating conditions — the two periods coexist (see figure 

It is noted that at higher bias voltages the number of peaks in the differential conductance 
increases, as new current paths become available at different energies for different occupation 
numbers. Only when the ratio U/A^ is an integer do peaks corresponding to different 
occupation numbers coincide. 

The Ohmic conductance through a quantum dot (figure ^) differs from that through a 
metallic dot in two significant ways. Firstly, the periodicity of the conductance peaks has 
increased by the level spacing A^. Secondly, the temperature dependence of the peaks has 
changed its nature. An increase in temperature now not only leads to larger thermal broad- 
ening, but also to a lowering of the peak amplitude which is inversely proportional to the 
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temperature. This is due to the fact that transport proceeds through a single energy level. 
The temperature dependence is proportional to the derivative of the Fermi-distribution func- 
tion in the reservoirs^. Therefore, at low temperatures k^T « U the Ohmic conductance 
can be written as 



^^^^ ~ hT, + TnV AksT cosh\^) ^^^^ 
where the successive charge degeneracy points are spaced by an energy U +A^. For higher 
temperatures several levels should be taken into consideration, each of which is weighted by 
a factor determined by the Boltzmann distributionS. 

where Pi{N) is the joint probability that the dot contains electrons and that the single 
particle level is empty. Note that the contributions from levels other than at positions /iat 
are too small to cause a peak in the conductance. 

Apart from the I-V characteristics and the Ohmic conductance, there is a third useful 
experiment that can be carried out. This involves applying a constant bias across the dot 
and studying the resultant current as a function of the gate potential. Typically the source- 
drain bias is chosen to be less than the charging energy which isolates the effects of the 
energy separation of the zero dimensional states of the dot. In figure ^ the current has 
been calculated as a function of the gate voltage for a dot with a constant level spacing 

= O.IU. The thermal energy is given by ksT = O.Olf/. In accordance with some recent 
experimental^ a number of peaks and troughs can be observed in the current. All peaks 
and troughs can be characterised by the number of levels available to an incoming electron 
to tunnel onto, and the number of levels from which an electron can tunnel out of the dot. 
For a Fermi level separation 6fi of 0.15U between the reservoirs (see figure |^) these numbers 
are 1, 2 — > 1, 1 — > 2, 1. This explains why the total peak is split into two subpeaks. For 
the second set of graphs with 6fi = 0.27U the sequence of available tunnelling levels is 
1, 3 ^ 1, 2 ^ 2, 2 -H> 2, 1 3, 1. It is also clear that asymmetric tunnelling barriers cause 
the current peaks to be asymmetric. 
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IV. TUNNELLING THROUGH TWO DOTS IN SERIES 



When studying the current and conductance properties of two quantum dots connected 
in series between the source and the drain, one obviously has to take into account the 
tunneUing between the dots. Inter-dot transitions are quahtatively different from transitions 
between a dot and a reservoir. The reservoir can be assumed to have a continuous density 
of states, so that electrons can always tunnel elastically into and out of the reservoir. Even 
though inelastic tunnelling events would in principle be allowed, their contribution would be 
relatively small compared to the elastic tunnelling rate. Therefore inelastic scattering only 
has a significant effect on the transport through a single dot if the scattering takes place 
inside the dot, i.e. relaxation. 

As far as tunnelling between dots is concerned, it is clear that the elastic tunnelling rate 
is significant only when the energy levels in the dots line up. This is obviously not generally 
the case. Usually an electron would have to interact inelastically in order to tunnel to a 
different energy level in the neighbouring dot. The energy difference would normally be 
absorbed or provided by phonons. 

As a theoretical model, one can consider the double dot system to be coupled to a 
phonon reservoir or a heat bath with a coupling strength The energy spectrum 

of the independent oscillators of the phonon reservoir is characterised by the density of 
phonon states pph{E). Some work has been done to calculate the current through a double- 
barrier resonant structure with some interaction between electrons and photons&H. The 
coupling to the optical phonons creates transmission subbands which are observable in the I- 
V characteristics. In this section, where inter-dot transmission is considered, the interaction 
with acoustic phonons will be the dominant mechanism. The energy spectrum for acoustic 
phonons is given by the Debye density of states. 



In a two-dimensional electron gas the Debye temperature 0d is approximately in the range 




(27) 



11 



Qd ~ 200 - 700kE. This is of the order of a few tens of meV, which corresponds to several 
times the Coulomb interaction energy U. This is typically larger than the voltage drop 
across the device, so that one can simply use the phonon density of states below the cut-off 
energy ksOD- 

Phonons are bosons, so that the Pauli principle does not apply and states can be multiply 
occupied. The average occupation number of states at a given energy is given by the Bose- 
Einstein distribution nsiE). 

e{E)-ei-E) 

where Q{E) is the Heaviside step function and E{—E) is the energy provided (absorbed) 
by the phonon bath. The Bose-Einstein distribution at positive and negative values of E 
differs by 1, which is indicative of the fact that a heat bath can always absorb energy. 

A plausible model for the inelastic tunnelling rate T^2 between an energy level in the 
first dot and ei in the second dot is as follows 



|2 



I , I ^ foo , ^ 



U, ,'m2 



Pph{uj' - uj)nB{uj' - uj) (29) 



where Pi{uj) and P2(^') densities of state for the energy levels. Considering the dots 

in isolation from the reservoirs, the eigenstates cannot be regarded as localised in each of 
the dots, but the wavefunctions will leak slightly into the other dot by virtue of the hopping 
potential Vm. This results in the overlap integrals \(f)'l{uj)(t>[{Lj')\'^ and |02(t^)02(^')P- Naively 
this can be interpreted as the probability that the inelastic process can take place between the 
wavefunction components in the same dot. The wavefunctions are calculated from a simple 
two-dimensional Hamiltonian with diagonal elements u and u' and off-diagonal elements 
given by the hopping potential Vm- By calculating its eigenvectors the overlap matrix can 
be determined. 

i:i/M^-(.-)r^ (,,_yV4va (^°' 

The density of states in the dots is broadened mainly as a result of the inelastic scattering, 
which is a prerequisite for the incoherent regime. To a first approximation the broadening 

12 



can 



be considered as a Lorentzian distributionii with a broadening proportional to the 



phase-breaking rate. Now equation 29 can be rewritten as a single integral over the energy 



difference E = uj' ~ uj between the initial and final state. 

- V /_: <'^(^T^^^^^".(^) (31) 

with Aph = 2Vlj\Vph\^Cph/Ti. 

It must be stressed that the broadening of the energy levels should be much less than 
the temperature. Otherwise it is not valid anymore to use the Boltzmann distribution to 
calculate the occupation probabilities. This assumption is consistent with the calculation of 
the reservoir-dot transition rate of section |T|. Figure ^a shows the energy dependence of the 
inter-dot tunnelling rate for the case when the energy-levels are considered to be (5-functions. 
However, it should still be taken into account that the level broadening is large compared 
to the coupling between the dots, causing the peaks in the transition rate at = ^2Vm to 
be smeared. At energy differences A = — which are small compared to the broadening, 
the transition rate can now be approximated by a Lorentzian of width 2T^/h and height 
2Ap/iln[r0/2Vjvf]/7rr(^ (see figure |^). At larger energy differences |A| > F the function is 
sufficiently close to the Bose-Einstein distribution. 

As the broadening is much smaller than the thermal energy, the inter-dot transition 
rate for small A can get much larger than the transition rate between dot and reservoir so 
that the inter-dot transition is no longer the current-limiting process. In this case it is of 
no relevance whether the transition rate between the dots is infinite or simply very large. 
This suggests that the inter-dot transition rate may be approximated by the Bose-Einstein 
distribution AphnB{^) at all values of A. This makes the mathematical analysis much more 
transparent. 

When one calculates the transition rate between dots with a given occupation number, 
one should sum over all possible tunnelling events according to equation ||. The greatest 
contribution should come from the energy range where the electron can tunnel from levels 
which are mostly occupied to levels which are mostly empty. Outside this window the 
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product Pi(A;|A^i + 1)[1 — P2{l\N2)] falls off exponentially. However, since the Bose-Einstein 
distribution diverges at A ^ 0, the total transition rate seems to be dominated by matching 
levels, even when they are situated at energies far removed from the Fermi level. This 
is clearly an unphysical situation. This anomaly can be removed by again including the 
broadening in the calculation of the transition rate. Alternatively, it can be argued that 
the transition rate for a given pair of levels is the combination of the previously defined 
rate and the rate at which electrons can get into an excited level. This second rate is the 
relaxation rate h/Tj-ei- Remembering that the inverses of the rates of consecutive processes 
add together, the total inter-dot rate can be obtained for given occupation numbers. 

TuiNu = E p^{m + m- p^±im.±i)] 3^'^T (^2) 

k,l Jl2 + Al/T-rcl 

In order to perform some realistic simulations it is helpful to know the dependence on the 
size of both the charging energy U and the confinement energy A^. In the charging model 
approximation the Coulomb repulsion is inversely proportional to the capacitance C of the 
dot. This implies that the repulsion energy U is also inversely proportional to the area of 
the dot. The single particle energy level spacing is given by0 

= — 33 

m area 

where m* is the effective mass of an electron in the two dimensional electron gas in which 
the quantum dot has been defined. It follows that both the Coulomb energy and the level 
spacing scale inversely proportionally to the area of the quantum dot. 



A. Ohmic conductance 

Figure |^ shows the Ohmic conductance through two dots with a relatively strong inelastic 
tunnelling coefficient Ap^. In order to analyse the structure of the peaks, assume that the 
total occupancy of each dot can only fluctuate by one electron. Moreover, assume that 
only a single level per dot (at the charge degeneracy points) contributes substantially to the 
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transport. Then the global master equation is used to obtain an expression for the Ohmic 
conductance. 



AksTcos\.{^) cosh(f^) 



r|coshf^) r^cosh 



cosh cosh ( 



+ 2AphfiVLVR cosh ( ^^r^) + 2VlTr{Vl + Tr) sinh ' ' 




(34) 



2kBT J V 2A;bT 

When the two reservoirs are equally strongly coupled to the dots, i.e. = = F, then 
the above expression peaks at /i = (e^ + e/)/2. The peak conductance is given by (using 
A = e/ - ek) 



SksT cosh^ (i^) cosh (4^) + 2F sinh 



|A| 



(35) 



When the peak height is investigated as a function of the temperature, then it appears that 
it has a maximum at a value of A;bT/|A| which is given by the solution of the following 
transcendental equation. 

V r keT kgT I yiksT) T keT ^ ' 

The solution is plotted in figure ^. When the ratio of the inelastic tunnelling coefficient to 
the reservoir coupling Aphh/T is of order unity or larger, then the temperature at which a 
particular conductance peak reaches its maximum height is given by roughly half the energy 
difference |A|. This maximum height is larger as the energy difference A gets smaller. 
This calculation is in good quantitative agreement with the temperature dependence of 
the conductance curves of figure 0. At higher temperatures the calculation becomes more 
inaccurate as several levels and more occupation numbers will start to contribute to the 
transport. 

In figure ^ the conductance is shown for a double dot system which is identical to that of 
figure with the exception that the inelastic tunnelling coefficient is two orders of magnitude 
smaller. From figure ^ one would expect the peaks to be maximised at a temperature 
ksT ~ 2.2|A|. The above description of the conductance peaks seems to apply to most 
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of the peaks. However, it is clear that the peak which is situated at /i ^ —4.711 has an 
anomalous behaviour. The conductance at this point is much larger than expected. This is 
due to the fact that one of the lower levels in the first dot very nearly matches up with the 
dominant level in the second dot, thus strongly enhancing the inter-dot tunnelling rate. 

This effect most strongly shows up when the inelastic tunnelling coefficient Ap^h is small 
compared to the coupling F to the reservoirs. This is simply due to the fact that the inter- 
dot tunnelling is the main current-limiting process and will therefore tend to dominate the 
physics. For large inelastic tunnelling coefficients the current will mainly be determined by 
the matching of the levels between dot and reservoir. This explains why the afore-mentioned 
effect is almost imperceptible in figure ^. 



B. Current characteristics 

The current-voltage characteristics have been calculated in figure ^ for a double dot system 
with specifications as indicated in the caption. Figure ^ shows the results for a system 
where the dot specifications have been swapped. Physically this simply amounts to varying 
the chemical potential in the right reservoir instead of the left reservoir. 

At low temperatures two periods may be observed in the I-V characteristics. The larger 
period is given by the Coulomb repulsion in the first dot plus a single particle spacing 
Ui + Ag^. The smaller period is simply given by the single particle spacing A^^. This 



behaviour is reminiscent of the current through a single dot with F^, ^ Tr (see section |T|). 
In other words, the second dot with its surrounding barriers acts as a single barrier with a 
reduced transparency. The electrons which have entered the first dot will tunnel into a lower 
energy level in the second dot at a rate of approximately Aph since the levels in the dots will 
not normally line up. At larger temperatures the current curves will lose some features due 
to thermal smearing. This is clearly the case in figure pT[ 



However, the current curve of figure 10 has a more complicated form at high tempera- 



tures. It has a region of negative differential conductance which occurs in the bias range 
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where the average occupancy of the first dot increases from one to two electrons with re- 
spect to the average occupation number at zero bias. This can be explained by considering 
the (unlikely) case where a pair of energy levels matches up for a given set of occupation 
numbers. This will result in a highly increased tunnelling rate between the dots. When 
the bias across the device is raised the average occupancy of the dot will increase. This 
means that the dots are less likely to contain the number of electrons for which the levels 
lined up. Consequently the tunnelling rate between the dots and hence the overall current 
will decrease, in spite of the fact that the total number of tunnelling paths is likely to have 
increased. 

The temperature at which negative differential conductance can occur (if at all) is set 
by the energy difference of the matching pair of levels in the dots. Their energy separation 
has to be significantly less than ksT for the tunnelling rate to increase dramatically. A 
higher temperature can cause some levels to match up which could not really be considered 
energetically close at lower temperatures. However, a higher temperature also means that 
there is a smaller probability that the total inter-dot tunnelling is dominated by just a 
single pair of levels. This will reduce the effect. The conclusion is that negative differential 
conductance is more likely to occur at higher temperatures, but when it occurs at a lower 
temperature the effect will be more pronounced. 

Similar to the one dot case, one can extract information about the level spacing when one 
considers the current through the system at fixed bias 6fi while varying the gate voltage of 
one of the dots, a measurement first performed by van der Vaart et a/.0 This is depicted in 



figure |T^. The Coulomb repulsion energies and the level spacings are the same as before. The 
inelastic tunnelling coefficient is small compared to the dot-reservoir coupling. The empty 
levels are the levels that an excess electron can occupy. Note that no more than a single 
extra electron can be contained in the dot, since this is prevented by the Coulomb blockade. 
Since electrons carry a negative charge, a rise in the gate potential causes a downward shift 
of the energy levels. 



The results are shown in figure |1^ for a range of values of Sjj. A few generic features 
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can be noted which are also true for tunneUing through a single dot. Firstly, the current is 
periodic in the gate voltage with a period Ui + A^^. Only a single period is shown in the 
figure. Secondly, a significant current is only allowed to fiow when the Coulomb blockade is 
overcome in both dots. This requires the lowest available level in each dot to lie within the 
energy window 6fi. This is clearly shown in figure [T^ where a larger bias voltage allows the 
current to fiow at more values of the gate potential. 

The most striking feature of figure ^ is the occurrence of sharp current peaks. These 
happen at values of the gate voltage where one of the occupied states of the first dot (con- 
taining the excess electron) lines up with an empty level in the second dot. It is clear that 
this will happen with a periodicity A^^ (in the figure A^^ = 0.13f/). When the bias voltage 
is such that several levels in the second dot are contained within the energy window 6^, 
then peaks will also occur at intervals A^^ ■ This is the case in the third graph of figure p!3| , 
where A,^ = 0.32f/. 

Finally, note that the current in the valleys between the peaks increases more or less 
linearly with the number of peaks. This refiects the fact that at a higher gate potential 
there are more electrons in the first dot which are able to tunnel downwards in energy into 
the second dot. Each tunnelling path has an associated off-resonance tunnelling rate of Ap^ 
which will contribute an amount Aphhio/T towards the total current (assuming Aphh <^ F). 
Therefore the current minima should increase by this amount every time the gate potential 
is increased by an amount A^^ . 

This seems to be a much more powerful method for determining the single particle level 
spacing than the analogous experiment with a single dot. In the model used in this chapter, 
the line shape will be a multiple of the Bose-Einstein distribution. However, at very small 
energy separations A the intrinsic width of the levels becomes significant and the lineshape 
will be approximately a Lorentzian of width 2F(^ near resonance but will be asymmetric 
off-resonance. This seems to corresponds at least qualitatively to recent experiment^. In 
the case of large Api^h/T, the middle barrier is no longer the current limiting barrier and 
the current will not peak anymore but simply increase with A until a saturation value of 
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^e^/h is reached. Further experiments will have to show whether it is justified to assume 
that the inelastic tunnelhng rate can be approximated by the Bose-Einstein distribution. If 
this turns out to be a false assumption, then the master equation can still be used to model 
the effect of a more realistic inelastic tunnelling rate. 

V. CONCLUSIONS 

For a dot with discrete levels, the distribution function at equilibrium differs from the Fermi- 
Dirac distribution, especially at low temperatures. The peaks in the Ohmic conductance 
are spaced by an amount U + A^. In the I-V characteristics the effect of both the charge 
quantisation and the size quantisation can be observed. 

The transport through a double dot has been investigated assuming that inelastic trans- 
port through the inter-dot barrier can take place by means of interaction with acoustic 
phonons. The main peaks in the Ohmic conductance reach a maximum at a temperature 
ksT which is at least roughly half the energy difference A between the dominant levels in the 
two dots. Other levels in the two dots that are well ahgned can also have a significant effect 
on the conductance, especially when the inter-dot coupling is weak. The I-V characteristics 
may contain regions of negative differential conductance. This is less likely to occur at low 
temperatures, although its effect will be stronger than at higher temperatures. 

The inter-dot spacing can be analysed spectroscopically by investigating the current 
through the double dot at fixed bias as a function of one of the gate voltages. This produces 
very narrow peaks with a width that is closely related to the intrinsic level width. This 
method eliminates thermal smearing of the peaks. 
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FIGURES 
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FIG. 1. Distribution function for a dot with equally spaced energy levels 




FIG. 2. I-V characteristics for various level spacings (F^ = Tr) 
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FIG. 3. Differential conductance for various level spacings {Tl = Tr) 



15 



10 - 



"D 
C 

o 
o 



kT = 0.01 U 
kT = 0.02 U 
kT = 0.05 U 



/ 



\ 



-1 1 

chemical potential (U) 



FIG. 4. Ohmic conductance through a dot with discrete levels (Ag = O.IU) 
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FIG. 5. Current vs gate voltage at fixed bias (A^ = O.IU) 
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FIG. 6. Energy dependence of the inter-dot transition rate: (a) ignoring level broadening, (b) 

including broadening = 5Vm- 
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FIG. 7. Ohmic conductance through a double dot with Aphh/T = I (Ui = OAIU, U2 = U, 
Agj = 0.13C/, = 0.32C/). The empty and filled circles indicate the positions at which the 
average occupation increases by one for dot 1 and 2 respectively. 




FIG. 8. Temperature at which a conductance peak reaches its maximum 
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FIG. 9. Ohmic conductance through a double dot with Aphh/T = 0.01 {Ui = OAIU, U2 = U, 
Agj = 0.13C/, Agj = 0.32[/). The empty and filled circles indicate the positions at which the 
average occupation increases by one for dot 1 and 2 respectively. 
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FIG. 10. I-V characteristics for a double dot at various temperatures {Ui = 0.41J7, U2 
= 0.13?7, Ae^ = 0.S2U, Aphh/F = 0.01) 
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FIG. 11. I-V characteristics for a double dot at various temperatures (the dots are interchanged 
with respect to the previous figure) 




FIG. 12. Energy diagram for tunnelling through two dots at fixed bias 
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FIG. 13. Current vs gate voltage of first dot at fixed bias {Ui = OAIU, U2 = U, A^^ = 0.13U, 
= 0.32C/, ksT = 0.02c/, = Tfl = F, Api^h/T = 0.01) 
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